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METHOD AND APPARATUS FOR DS-CDMA INTERFERENCE SUPPRESSION USING 

CODE-SPECIFIC COMBINING 

BACKGROUND OF THE INVENTION 

5 

The present invention relates generally to RAKE receivers and particularly to 
interference suppressing RAKE receivers. 

RAKE receivers represent a well-known approach to multi-path reception, 
particularly in Direct Sequence Code Division Multiple Access (DS-CDMA) wireless 

1 0 communication. In a multi-path communication channel, a transmitted signal travels through 
multiple propagation paths to the receiver. Thus, the receiver receives multiple "echoes" of 
the transmitted signal, with each multi-path echo generally suffering from different path 
delay, phase, and attenuation effects. 

RAKE receivers exploit multi-path propagation by assigning each of two or more 

1 5 RAKE "fingers" to one of the incoming multi-path echoes. Each finger is tuned to a 

particular multi-path echo. By estimating the channel effects, e.g., phase and attenuation, and 
by properly accounting for the differences in path delays, the individual output from each 
RAKE finger may be combined with the outputs from the other fingers to provide a 
combined RAKE output signal with a greatly improved signal-to-noise ratio (SNR). 

20 In DS-CDMA systems, such as Wideband CDMA and IS-2000, high transmission 

data rates are achieved by transmitting data at a low spreading factor and/or on more than one 
spreading code (multi-code). When a low spreading factor and/or multi-code is used, 
performance is sensitive to multi-path dispersion. With dispersion, there are multiple echoes 
of the transmitted signal with different relative delays. These echoes interfere with one 

25 another. Not only is orthogonality lost between successive symbols as one symbol overlaps 
with the next, but orthogonality is also lost between symbols sent on different, orthogonal 
codes. 
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As a result, performance is often limited by interference between different symbols 
being sent to a particular user. These symbols can correspond to the previous or next symbol, 
symbols sent in parallel on another code carrier, or both. In general, this interference is 
referred to as self-interference, where self-interference may include intersymbol interference 
5 (ISI) or inter-code interference (ICI). A key aspect of ISI is that it varies from symbol to 

symbol. This variation is due to the spreading codes being a combination of Walsh codes and 
a common scrambling code. The scrambling code has a period much longer than the symbol 
period, which effectively changes the overall spreading code from symbol period to symbol 
period. 

10 

SUMMARY OF THE INVENTION 
The present application discloses a method and apparatus for reducing interference, 
such as inter-symbol interference, from a symbol of interest in a RAKE receiver. The RAKE 
receiver comprises a plurality of RAKE fingers, a processor, and a combiner. The plurality 
15 of RAKE fingers despread symbols received over multiple paths of a multi-path channel. 
The processor determines cross-correlations between different symbols. The combiner 
combines the despread symbols using the cross-correlations to reduce interference from the 
symbol of interest. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 comprises a block diagram of a conventional generalized RAKE (G-RAKE) 
receiver. 

Figure 2 comprises a block diagram of an exemplary multi-symbol RAKE (M-RAKE) 
receiver according to the present invention. 
25 Figure 3 comprises a block diagram of an exemplary restricted M-RAKE receiver 

according to the present invention. 
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Figure 4 illustrates an exemplary multi-channel linear transversal filter for the 
restricted M-RAKE of Figure 3. 

Figure 5 illustrates an exemplary linear transversal filter for the multi-channel linear 
transversal filter of Figure 4. 

Figure 6 illustrates an exemplary linear equalizer according to the present invention. 

Figure 7 illustrates an exemplary multi-code linear equalizer according to the present 
invention. 

Figure 8 illustrates an exemplary multi-channel filter for the multi-channel linear 
equalizer of Figure 7. 

Figures 9A and 9B illustrate a mathematical technique for simplifying weight 
calculations associated with RAKE receivers. 

Figure 10 illustrates an exemplary wireless communication system that utilizes the 
RAKE receivers of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
Figure 1 illustrates a conventional wireless communication receiver 100 utilizing a 
generalized-RAKE (G-RAKE) 112. Wireless communication receiver 100 comprises 
receiver front end 1 10 and G-RAKE 1 12. The receiver 100 and/or the G-RAKE 1 12 may be 
embodied in an application specific integration circuit (ASIC). Receiver front end 110 
outputs received signals r(t) to G-RAKE 112. These received signals r(t) comprise streams 
of sample values obtained from wireless signals received by one or more antennas associated 
with one or more receiver front ends 110. An exemplary front end 110 may include 
amplifiers, filters, mixers, digitizers, and/or other electronics as needed to produce a sampled 
signal suitable for processing by the G-RAKE 112. The received signal r(t) output to G- 
RAKE 1 12 typically comprises a composite signal that includes one or more multi-path 
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echoes arising from multi-path propagation and which arrive at the receiver 100 from 
different directions and with different time delays. Further, received signal r(t) may include 
components transmitted from or received by different antennas (not shown). The task of the 
G-RAKE 1 12 is to generate an estimate s(m) of each transmitted symbol s(m) from the 
5 received signal r(t). 

G-RAKE 1 12 comprises a plurality of RAKE fingers 120, RAKE processor 130, and 
RAKE combiner 132. Each RAKE finger 120 processes different time shifts or multi-path 
echoes of the received signal r(t). Typically, each RAKE finger 120 comprises a delay 
element 122 and a correlator 124. Delay elements 122 delay the received signal r(t) to time 

1 0 align the multi-path echoes processed by each RAKE finger 120. Correlators 124 correlate 
the delayed signals with a spreading code to extract the assigned multi-path echoes from the 
received signal r(t), Despread values from correlators 124 are combined in combiner 132. 
Combiner 132 typically includes weighting elements 126 and summer 134. Weighting 
elements 126 weight the multi-path echoes output from respective correlators 124. The 

1 5 weighted multi-path echoes are summed symbol-by-symbol by summer 1 34 to form an 

estimate of a symbol of interest, s(m) 9 during each symbol period. While some of the figures 
refer to the estimated symbol of interest as s(m) 9 throughout this specification the symbol of 
interest s(m) and the estimated symbol of interest s(m) may be represented by s(0) and §(0), 
respectively, where w=0, without loss of generality. Further, those skilled in the art will 

20 appreciate that the time shifts associated with delay elements 122 correspond to the delays of 
the multi-path echoes, while the combining weights associated with weighting elements 126 
may correspond to the conjugates of the multiplying coefficients of the multi-path echoes 
(RAKE) or conjugates of weights that depend on the coefficients and a noise correlation 
matrix (G-RAKE). 
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In the conventional receiver 100 described above, the symbol estimates s(m) output 
during each symbol period may be corrupted by interference, such as intersymbol 
interference (ISI), which may result in demodulation and decoding errors. The ISI may be 
due in part to the dispersive effects of the communication channel, and in part to cross- 
5 correlations between the spreading codes used for different symbols. It should be noted that 
the spreading code is a combination of an assigned Walsh code and a scrambling code. The 
period of the scrambling code is typically larger than a symbol period, so the net effect is a 
time-varying spreading code. The present invention uses knowledge of the communication 
channel and cross-correlations between spreading codes for different symbols to determine 

1 0 cross-correlations between symbols to reduce interference. 

Figure 2 illustrates an exemplary multi-symbol RAKE (M-RAKE) 200 according to 
the present invention. As with the G-RAKE 112, M-RAKE 200 may be embodied in an 
ASIC. M-RAKE 200 comprises a plurality of RAKE fingers 220, a RAKE processor 230 
and a RAKE combiner 232. Unlike the RAKE fingers 120 of G-RAKE 1 12, RAKE fingers 

1 5 220 may despread symbols from different symbol periods. Typically, one or more RAKE 
fingers 220 will have a delay that corresponds to a symbol of interest s(0). Other RAKE 
fingers 220 will have delays that correspond to interfering symbols (i.e. symbols near the 
symbol of interest that contribute to intersymbol interference). During each symbol period, 
weighting elements 226 weight the despread symbols output from corresponding correlators 

20 224 to form weighted outputs, which are then summed in summer 234 to form a symbol 

estimate s(0). RAKE processor 230 calculates the combining weights of weighting elements 
226 to cancel or reduce intersymbol interference due to the effects of the communication 
channel and cross-correlations between spreading codes for different symbols. 

The combining weights or weighting factors for M-RAKE 200 may be determined 

25 from a maximum likelihood (ML) or a minimum mean square error (MMSE) solution 
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computed by RAKE processor 230 utilizing knowledge of the different symbol spreading 
sequences and the channel coefficients. Typically, these combining weights are based on an 
analysis of the correlator outputs. With many of the mathematical details and definitions 
deferred till later, the vector of correlator outputs may be represented as: 

5 y= [y( rf i>™i) y( d i> m 2) - y( d j>™j)Y 

=h 0 *(0) + £iv(o + n, (i) 

where y(d p corresponds to the output of delay element J and correlator J, 

h 0 s(0) + y^h ; ^(i) represents the vector output components associated with each transmitted 

1*0 

symbol, and n represents a noise output vector (thermal noise, other interference). Generally 
1 0 speaking h, comprises the composite channel impulse response associated with the i th symbol, 
s(i) 9 which includes the channel effects and despreading (spreading sequence cross- 
correlations). The first term in Eq. (1) represents the symbol of interest and the summation in 
Equation 1 represents the contribution of all other symbols, or the ISI, in y. 

The symbol of interest, s(0), is estimated by using a linear combination of the outputs 
1 5 of RAKE fingers 220, which may be represented as s(0) = w H y , where is the conjugate 
transpose of a weight solution vector w. The linear combination is performed in RAKE 
combiner 232 by weighting elements 226 and summer 234. 

The weight solution vector, w, based on a maximum likelihood criterion for the 
estimate, may be calculated by: 

20 w = R u h 0 , (2) 

where R M is the correlation matrix of all undesired components in Equation 1 . Correlation 
matrix R u may be calculated by taking the expectation with respect to noise and random 
transmitted symbols (assuming the random transmitted symbols are independent). The 

25 resulting correlation matrix of the undesired components may be represented by: 
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R M =EMr+R« 5 (3) 

where R„ is the correlation matrix of the noise, which may be correlated across fingers after 
the despreading process. The calculation of the noise correlation matrix R„ is shown below. 
Because the expectation is taken over random transmitted symbols, instead of the known 
5 spreading codes, the combining weights applied by weighting elements 226 will change each 
symbol period due to the time- varying nature of the spreading waveform. 

The derivation of the correlator outputs defined in Equation 1 will now be described. 
The transmitted signal for a single code channel may be defined as: 

x(t) = jEf d s(i)a t {t-iT) i (4) 
1 0 where E is the symbol energy, and s(i) and aft) are the i th symbol and its spreading 

waveform, respectively. The spreading waveform aft) is composed of a complex sequence 

of code chips convolved with a pulse shape. It is assumed that the spreading waveform is 

normalized so that J° |a i (/)| 2 *fc = 1 • 

If the multi-path channel is discrete and comprises L rays, then 

15 g(0 = i d g,S(t-r l ) (5) 

represents its impulse response where gi and r, are the complex gain coefficient and delay of 

the I th ray, respectively. The signal received through this channel may be represented by: 



r(/) = V^Z Z it - iT - t,) + ff (0 , (6) 
where n(t) is a complex noise process. 

20 The finger of a RAKE receiver incorporates a delay, d, which precedes a correlator 

matched to the spreading waveform of the m th symbol, a m (t). The output of the finger may be 

given as: 

y(d,m)= f r(a + d)a m (a)da 

J— CO 
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L oo 

= ^ L Z Z {a + d-iT- r, )a m (a)da + « (rf, m) , (7) 
where n(d,m) is the complex noise output of the correlator. Assuming that the symbol of 

interest is 5(0), the traditional RAKE would have m = 0 and the delays of the fingers would 

be matched to the delays of the resolvable rays of the multi-path. However, for the M-RAKE 

200 of the present invention, the spreading waveform used for correlation may be generalized 

to be any symbol's spreading waveform. Equation 7 may be rewritten as: 

y(d,m) = £j(a + d)a m {a)da 

= ^ L Z Z {a + d-iT- r t )a m (a)da + n (d, m) 



10 = J /*,(</ -iT,m)s(i) + n(d,m) (8) 

where 

« m) = f «(a + t)a * (a)rfa (9) 

is the colored noise output of the correlator, and 

R ^ = £fl|(a + r)a;(a)i/a (10) 

15 is the correlation of the i th symbols spreading waveform against the m xh symbol's spreading 
waveform, used by the finger's correlator. Finally, the components associated with a 
particular symbol (a correlation between symbol m and i) may be grouped together into 

W>™) = lLgA, m (t-T l ) 9 (11) 

a net channel impulse response for the i h symbol, where spreading code and pulse shape 
20 information provides Ri >m (t) and channel estimation is used to get g/. 
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Finally, the RAKE finger outputs may be vectorized for the detection of the symbol of 



interest as: 



y = [y(«f. 



(12) 



where there are / RAKE fingers 220. In a traditional matched filter case, J = L, the number 
5 of resolvable rays in the multi-path channel. Also, J> L fingers may be used. 

The calculation of the noise covariance or correlation matrix R„ referenced in 
Equation 3 is computed in one of two ways. In the first way, the noise correlation matrix 
represents white noise, n(t) 9 which has been colored by the correlation process (filtering) of 
the RAKE fingers 220 to produce n(d,m), where the noise is a function of the delay d of the 
1 0 RAKE finger 220 and the particular spreading waveform of symbol m used by the correlator 



224. 



Defining the noise output vector of the RAKE fingers 220 as 



n= [n(d X9 m x ) n(d 29 m 2 ) 



n(dj,mj)] T 9 



(13) 



the noise correlation matrix may be represented by: 



15 




(14) 



The cross-correlation of the ( finger and the f finger may be given as: 



R„ (/, j) = E${di . m i )» * ( d j > m j)} 





(15) 



20 Because the noise is white, E{n(t)n *(t + r)} = a 2 8(t - r) , where a 2 is the noise power, 



Equation 15 becomes: 



R - = L L < < a K iP^{a + d,-0- dj)dad/3 . 



(16) 



Sifting property of the S( ) function in integration results in: 
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R n 0'»y)=o- 2 £<*o(«K(«+4 ~dj)da 

-cr'Rm^-dj). (17) 

Thus, the noise covariance can be constructed using spreading code cross-correlation 
information and noise power estimates (a 2 ). The latter can be estimated by using known 
techniques. When interference from other user signals is present (i.e., not modeled as ISI or 
ICI), then the noise at the receive antenna may not be white. A second way to compute the 
noise covariance matrix is to estimate noise correlations using despread values. For example, 
the model in Equation 1 also applies to the pilot symbols being sent on the pilot channel. 
This model can be used to obtain vectors of noise samples, n, by subtracting the pilot symbol 
contributions from the y vector. Outer products of the noise sample vectors can be formed 
and smoothed in time to give an estimate of the noise covariance matrix. As this matrix is 
Hermitian symmetric, only the upper or lower triangle of the matrix needs to be estimated. 

Turning now to Figure 3, an alternate embodiment of the present invention, referred 
to herein as a restricted M-RAKE 300, will now be discussed. Restricted M-RAKE 300 
comprises a plurality of RAKE fingers 320 and a multi-channel linear transversal filter (MC- 
LTF) 340. In restricted M-RAKE 300, each RAKE finger 320 includes a delay element 322 
and a correlator 324. In this embodiment, the delays associated with each RAKE finger 320 
are chosen so that the despread symbols output from each correlator 324 are time aligned. 
The despread symbols output from RAKE fingers 320 are input to the MC-LTF 340. As 
shown in Figure 4, the MC-LTF 340 comprises a plurality of parallel linear transversal filters 
(LTFs) 342, each associated with a corresponding RAKE finger 320, and a summer 344 that 
adds the filtered outputs of each LTF 342 to generate a symbol estimate s(m) . Those skilled 
in the art will appreciate that other embodiments of the MC-LTF 340 are possible. In 
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general, past finger outputs are stored and results from multiple fingers and multiple symbol 
periods are combined. 

Each LTF 342 combines despread symbols output from a corresponding RAKE finger 
320 during a plurality of successive symbol periods as shown in Figure 5. The despread 
5 symbols combined by each LTF 342 include a symbol of interest and the nearby interfering 
symbols (both preceding and succeeding the symbol of interest) that contribute to ISI. The 
function of the LTFs 342 is to reduce the ISI in the symbol of interest as described in greater 
detail below. The filtered signals output from each LTF 342, referred to herein as the filtered 
output signals, are added by multi-path summer 344 to generate the symbol estimates. The 

1 0 summer 344 performs a similar function as a summer in a conventional RAKE receiver to 
improve the signal to noise ratio of the final symbol estimates. 

Each LTF 342 comprises 2M+1 weighting elements 346, 2M delay elements 348, and 
LTF summer 350, where M represents the length of the ISI in symbol periods. Thus, when M 
equals 2, there will be four interfering symbols (two preceding and two succeeding) 

1 5 surrounding a symbol of interest s(m). LTF 342 receives and delays the despread symbols 
output from a corresponding RAKE finger 320 in a tapped delay line comprising a plurality 
of delay elements 348. The despread symbols may be corrupted by ISI. During each symbol 
period, the LTF 342 weights and combines the despread symbols output during 2M + 1 
symbol periods centered on the symbol of interest s(m) in weighting elements 346 and LTF 

20 summer 350. The weighting factor for each weighting element 346 is solved for jointly in a 
single step by processor 330 based on the channel coefficients and the cross-correlation 
between spreading codes, as described further below. 

Restricted M-RAKE 300 reduces the number of RAKE fingers 320 while providing 
the same RAKE finger outputs typically provided by a larger number of RAKE fingers 220 in 

25 the M-RAKE 200 shown in Fig. 2. For example, assume a restricted M-RAKE 300 
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comprises four RAKE fingers 320, where each RAKE finger output is fed into a multi- 
channel LTF 340 spanning ±1 symbols (M = 1) around the desired symbol. This exemplary 
M-RAKE 300 is equivalent to a non-restricted M-RAKE 200 with 12 RAKE fingers 220. 
Thus processing advantages may be accomplished with restricted M-RAKE 300. 
5 Figure 6 shows a linear equalizer (LEQ) 400 constructed in accordance with the 

present invention. LEQ 400 comprises a generalized RAKE (G-RAKE) 410 and a linear 
transversal filter 442 to combine successive outputs from the G-RAKE 410. Those skilled in 
the art will understand that the G-RAKE 410 may use a variety of combining methods, 
including traditional RAKE combining. G-RAKE 410 comprises a plurality of RAKE 

1 0 fingers 420, a RAKE processor 430, and a RAKE combiner 432. Each RAKE finger 420 
comprises a delay element 422 and a correlator 424. Like the restricted M-RAKE 300 of 
Figure 3, the delays associated with each RAKE finger 420 are chosen so that the despread 
symbols output from each correlator 424 are time aligned. The outputs of correlators 424 are 
weighted in weighting elements 426 and summed by summer 434 to form a RAKE output 

1 5 signal. The RAKE output signal from combiner 432 is input to LTF 442. As with the 

restricted M-RAKE 300 discussed above, the delay length of the LTF 442 is 2M + 1 symbols 
centered about the symbol of interest s(m). During each symbol period, LTF 442 combines 
2M 4- 1 RAKE output signals to generate a symbol estimate. The RAKE output signals are 
combined using weighting factors determined based on channel coefficients and cross- 

20 correlations between spreading codes to reduce ISI from the symbol of interest s(m). 

The weighting factors for the LTF 442 are determined according to essentially the 
same ML criteria way as with the M-RAKE 200, with some redefining of the data vectors in 
question. As the RAKE finger outputs are stored in LTF 442, the vector of RAKE finger 
outputs are redefined according to the time they were produced. That is, the outputs of 

25 correlators 424, previously defined in Equation 1, may be rewritten as: 
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y(mT) = \y(mT + d xy m) y{mT^d li m) ... y(mT + d J9 mj\ 

= h o (mrM0) + X h .( wr M0 + n(mr) s (18) 

where h/(ra7) is the response from the f h symbol seen in the m th RAKE finger output vector, 
and where the symbol of interest is represented by s(0). The output of combiner 432 
5 comprises the w th -combined output of the RAKE finger outputs and may be defined as: 

z(mT) = g H y(mT), (19) 
where, for example, the weighting and combining process simply comprises the complex 
conjugate of the coefficients of the channel taps in the dispersive channel, (traditional RAKE) 
g = [gi S2 gi] T - The contents of the LTF 442 may be represented in vector form as: 
10 z= [z(-MT) z(-(M-\)T) ... z(MT)f 

= h^(0) + Xh;*(0 + n', (20) 

Equation 20 has the same general form as Equation 1 , with the components of the various 
symbols and noise separated. The estimate s(m) of the symbol of interest may be 
represented as: 

15 i(0) = v"z. (21) 

The weight vector, v, in LTF 242 may be calculated by: 

v = R w " 1 h 0 (22) 
where the cross-correlation matrix R M . of the undesired portion of the signal may be given as: 

R u =XM: W +R„ (23) 

20 The cross-correlation matrix R M , describes cross-correlations between the effective 

spreading codes for different symbols. Details regarding the calculation of the noise 
correlation matrix R n , are discussed further below. As with the M-RAKE 200 of Figure 2, 
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the weighting factors calculated according to the ML criteria change each symbol period due 
to the change in the spreading waveform. 

Substituting Equation 18 into Equation 19 produces: 
z(mT) = g"y 

5 =g H h o (mT)s(0) + X g"^ (mT)s(i) + g H u(mT) , (24) 

from which the vector z may be written as: 



(25) 



15 





~ z{MT) ' 




' g w h 0 (MT) " 




" g w h,(Mr)" 




" g H n(A/r) " 


z = 








*((>) = £ 










z(-MT) 




g"h 0 (-MT) 


1*0 


g"h, (-MT)_ 







From this h' may be defined as 

h;=[ g H h,(Mr) g H h,((M-i)r) ... g»h ,(-mt)] t , 

1 0 and the noise vector, n', as: 

n' = [g H n(MT) g"n{{M-\)T) ... g"n{-MT)] . 
This noise vector n' may be rewritten as: 



n = 



g" 0" ■■■ 0" 

0" g H ... o" 



0" 0" 



g 



' n(MT) 
n(-MT) 



= G"n", 



where this block diagonal matrix of channel coefficients may be defined as: 



G = 



g 0 0 

0 g ••• 0 

0 0 g 



(26) 



(27) 



(28) 



(29) 
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This allows the correlation matrix R u . of the noise vector, n', to be written as 

= G H E{n n n nH ]G (30) 
= G H R n .G, 

Equation 30 shows that R„»may be calculated element by element in the same manner in 
which R„ was calculated (e.g., Equation 17 for the white noise case). 

The present invention may also be utilized in multi-code receivers. Multi-code 
transmission results when multiple code channels are devoted to a single downlink receiver to 
increase data carrying capability. This is analogous to several separate simultaneous 
transmissions being received by a single receiver. As the receiver is bound to demodulate 
each of these code channels, it can use the results of the demodulation of one code channel to 
aid in the demodulation of another code channel. Further, because of the multi-codes 
received by the multi-code receiver, each demodulated symbol may, in addition to ISI, be 
susceptible to inter-code interference (ICI) caused by symbols from adjacent codes. 

For example, the M-RAKE 200 of Figure 2 may be adapted to process multi-code 
signals by assigning a subgroup of RAKE fingers 220 to each code channel. While those 
skilled in the art will understand that tuning to the spreading codes of other code channels 
does not necessarily demodulate the other code channels, they will also appreciate that this 
tuning process provides information on the other code channels, which may help to reduce 
cross interference effects on the code channel of interest. In its most general form, the 
modified M-RAKE includes multiple subgroups of RAKE fingers 220 tuned to the spreading 
codes of different code channels. Further, these RAKE fingers 220 may be tuned to different 
symbol intervals and delays. 

The restricted M-RAKE 300 of Figure 3 may also be extended to multi-code by 
adding additional RAKE fingers 320 tuned to different code channels, with expansion of the 
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MC-LTF 340 according to the number of added RAKE fingers 320. In this scenario, the joint 
solution of the weights within the MC-LTF 340 takes into account the spreading codes of all 
code channels. 

Similarly, the LEQ 400 of Figure 6 may also be extended to multi-code by using a 
5 plurality of LEQs 400 tuned to different code channels and expanding the LTF 442 following 
each LEQ 400 to an MCO-LTF as shown in Figure 7. The multi-code linear equalizer 
(MCO-LEQ) 500 comprises three G-RAKEs 510, each tuned to one of three code channels, a 
RAKE processor 530, and a multi-code linear transversal filter 540. Each G-RAKE 510 
comprises a plurality of RAKE fingers 520 and a RAKE combiner 532. Each RAKE finger 

1 0 includes a delay element 522 and a correlator 524. The delays for the delay element 522 are 
chosen such that the output from each RAKE finger 520 is time aligned with the other RAKE 
fingers 520. Correlators 524 despread the received symbols, which are then weighted in 
weighting elements 526 and summed by summer 534. The weights applied by the weighting 
elements 526 are calculated in a conventional RAKE or G-RAKE manner as is well known in 

1 5 the art. The MCO-LTF 540 combines the RAKE output symbols from each G-RAKE 510 
using weighting factors determined based on cross-correlations between symbols based on 
channel coefficients and cross-correlations between spreading codes as previously described. 
In the situation where all of the code channels are allocated to a single user, the MCO-LTF 
540 outputs one symbol estimate for each code channel (co-channel). Thus, if there are three 

20 code channels, the MCO-LTF 540 would output three symbols per symbol period. In this 
case, the MCO-LFT 540 would reduce intersymbol interference attributable to symbols 
transmitted on the same code channel, as well as intercode interference due to symbols 
transmitted over different co-channels. The present invention, however, may also be used in 
the situation where each co-channel is allocated to a different user. 
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Figure 8 illustrates an exemplary MCO-LTF 540 corresponding to MC-LEQ 500 of 
Figure 7. The MCO-LTF 540 may, in one embodiment of the invention, comprise a plurality 
of multi-channel linear transversal filters (MC-LTF) 340 as shown in Figure 4. Each MC- 
LTF 340 combines RAKE output symbols from all of the G-RAKEs 510 using weighting 
5 factors determined based on channel coefficients and cross-correlations between symbols. 

Each MC-LTF 340 in the MCO-LTF 540 outputs a single symbol estimate corresponding to a 
selected one of the co-channels for each symbol period. The weighting factors for each LTF 
342 would differ. For example, the weighting factors for the MC-LTF 340 associated with 
co-channel 1 would be selected to reduce interference from the symbols on co-channels 2 and 
10 3. The same approach would be used for the MC-LTF 340 for co-channels 2 and 3, resulting 
in different weighting factors. Note that the LTF units that store the same RAKE outputs 
may share memory to reduce memory requirements. 

As discussed above, the weighting factors for LTF 342 may be determined according 
to an ML solution. However, computational savings may be achieved if the ML solution is 
1 5 calculated as a minimum mean squared error (MMSE) solution, and then converted to a ML 
solution through a scaling of the combining weights (or final output), which can be calculated 
from elements of the MMSE solution. 

It is straightforward to show that for a single user detection RAKE receiver 200, 300, 
400 the weight vector, w, ML combining of the RAKE finger outputs is the same as the 
20 MMSE weight vector solution, v, within a scaling factor. That is, 

w = av (31) 
where a is some real coefficient. This a coefficient essentially relates the reliability of the 
symbol estimate. 

Examining first the single user case, the RAKE finger outputs may be written as: 
25 y = h 0 j(0) + XM(0 + n, (32) 

i*0 
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where h,- is the net response across the fingers of the i th symbol, s(i), through spreading, 
channel, and despreading. The noise vector n represents the vectorized output of multiple 
fingers with white noise as the input. Without loss of generality, hos(0) corresponds to the 
symbol of interest where /w=0. 
5 For the MMSE solution, the data correlation matrix, = is{yy " }, may be used to 

solve for the familiar expression 

v = R; 1 h 0 , (33) 
where h 0 results from E{ys * (0)}, with this expectation taken only over the symbols. 

For the ML solution, the disturbance correlation matrix, R u , is used, which is the data 
1 0 correlation matrix less the part related to the desired symbol, or 

R M =R y -hX, (34) 
to solve for another similar, yet familiar expression 

w = R; 1 h 0 . (35) 

Equation 31 may be proven easily, starting with Equation 33 and using the fact of Equation 
15 34: 

v = R;'h 0 

= R ll l R B R > l lio 

= R;'[R,-Mo"fcX 

= k'hj-^hoj.h-^hj 

20 = w-wh^v. (36) 

Because v is a scalar, w may be written as: 
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resulting in a scaling factor of: 



1 



(38) 



a = 




This scaling factor, a, is written in terms of an MMSE solution, with the assumption that it is 
desirable to convert an MMSE solution to an ML solution. That is, the terms for calculating 
a would already have been be calculated for the MMSE solution. 

Other expressions for "a" can be used which account for the fact that the channel 
estimate may be a scaled version of the true channel coefficient, due to differences between 
the pilot and traffic channel power levels and spreading factors. Specifically, the factor "a" 
can be computed as 

a = — 

where A corresponds to the ratio of the energy in a traffic symbol to the energy in the pilot 
symbol. For example, A may be expressed as 



where rl is the ratio of power allocated to the traffic channel to the power allocated to the 
pilot channel. This ratio can be estimated or set to some nominal value. Estimation of this 
ratio is addressed in Bottomley et al., "Communications methods, apparatus, and computer 
program products using gain multipliers," U.S. Patent Application No. 09/968,443, filed 1 
October 2001 . Variable r2 represents the ratio of the spreading factor of the traffic channel to 
the spreading factor of the pilot channel (length of the pilot symbol used for channel 
estimation). These quantities are known at the receiver. 

Regarding the multi-user/multi-code detection case, which may be used in the 
WCDMA multi-code transmission, parallel data channels sent to one user are considered as 
multiple user signals. Separate banks of RAKE fingers, each tuned to a user spreading 
sequence, can produce a vector of RAKE finger outputs, z, and these may be combined 



A = tIt2, 



(40) 
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together to a vector symbol estimate, where the elements of the vector are the estimates of the 
individual user symbols. For the MMSE case, this may be represented by s MMSE = V"z and 

for the ML case, s MLSE = W"z where V and W are the MMSE and ML weight matrices, 

respectively. The same scaled relationship between the MMSE multi-user weight solution 
and the ML multi-user weight solution is not as obvious as in the single user case. 

Now consider the multi-code/multi-user case, where there are K codes in parallel. 
The RAKE finger output vector may be given by: 



K 



(41) 



j>0 A=l 
th % 2 th , 



where h,-,* is the net response of the k? user's t symbol, s(i, k). Again, n is the noise vector 
1 0 output from the fingers. 

Solving jointly for a vector of user symbols 

■(0 = [*(M) *0\2) 

produces: 



(42) 



1 MMSE 



(0) = 



[v, v 2 



v*r«, (43) 
where v* is weight vector to provide the MMSE estimate of the kf h user's symbol. This may 
be rewritten as a matrix V times the RAKE finger outputs 

W(0) = V"z, (44) 
where V = [vi v 2 . . . v K ] and s MMSE (0) represents the symbol of interest. It's straightforward 

20 to show that the joint MMSE solution that minimizes ||s(0) - s(0)|| 2 is the same as the solution 
that minimizes the squared error of the individual estimates of the user's symbols, 
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|s(0,/:) - £(0, k)\ 2 . As a result, the weight matrix for the joint MMSE solution, V, may be 
written as: 

v = [R;'h 0 , RX 2 ... R;-hJ 
= lR;'h w h 02 ... h 0 J 

= R^Ho, (45) 

where the net responses of the K users become the K columns of the matrix Ho, the 0 
subscript denotes that this is the response for the desired 0 th symbol vector. Going forward, 
the 0 may be dropped for clarity, understanding that the response vectors for the 0 th symbol 
are being considered. 

All of this vector and matrix manipulation provides a concise way of writing a joint 
MMSE solution. Assuming that this equation is solvable using matrix inversion, the 
implication of Equation 45 is that only a single matrix inversion needs to be calculated. For 
the ML solution, it is not obvious that this is the case. The vector ML solution may be 
written in much the same way as before 

s^(0) = W"z (46) 
where W = [wi W2 . . . Wk]. As a solution that gives us log likelihood ratios (LLR) for the bits 
of the individual users is desired, the estimates of the symbols to be ML individually is 
desired. 

The weight matrix, W, for the ML solution may therefore be represented by: 

w = [r>, r; 2 h 2 ... r;>J, (47) 

where the disturbance correlation matrix for each user's ML weight vector solution may be 
different, as denoted Ru * for the k th user. Keeping the assumption that the weights using 
matrix inversion are being solved for, it would appear that it is necessary to calculate K 
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different inverses to solve this system, as opposed to the single inverse matrix of the MMSE 
solution. 

However, it should be noted that as in the single user case, the disturbance correlation 
matrix for the ML solution is the same as the MMSE data correlation matrix less the 
component from the desired symbol, that is 

R uk =R,-hX. (48) 
It is recognized that the inverse of R W) * can be written as a "rank one update" inverse so, 



The h! h column of W becomes: 

10 w 4 = R: , h t + 



y k k y k 



= 

V 



h"v ^ 



(50) 



and the same conclusion as for the single user case is achieved, although shown from a 
1 5 slightly different direction (using the rank one update). Then Equation 47 becomes: 

W = [A,v, A 2 y 2 ... X K y K ] (51) 

where 

and ultimately where: 
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W = VA (53) 
where A is a diagonal matrix with X k as the kf h diagonal element. 

Grouping these individual solutions in a matrix equation provides a concise 
expression in Equation 53, but it also points out important implementation points. 
5 Specifically, given a solution involving matrix inversion, it emphasizes that only a single 
inverse matrix is calculated for an MMSE solution as opposed to K inverse matrices for the 
ML solution, which may represent significant processing savings. 

As discussed above and shown in Equation 54, the form of the MMSE solution is 
similar to the form of the ML solution. 

1 ° v^^Rj-'K (54) 

However, the correlation matrix of the data, R 2 is used, rather than the correlation matrix of 
the disturbance or noise, R M . 

The form of the z vector of RAKE outputs for the q th symbol may be represented by: 

z q =[z{q-M) z{q-{M-\)) ... z(g + M)f , (55) 

1 5 where the symbol period, T, is dropped for convenience. This vector essentially represents 

the contents of the LTF 442 corresponding to the LEQ 400 discussed above. In another form, 
z g may be rewritten as: 

z q =[z(q-M) aj r , (56) 

where most of the z vector is now represented by subvector a g . Note that for the next symbol, 
20 q + 1, the data is shifted through LTF 442, yielding: 

Vi =k zte + M + l)Y 9 (57) 
where z(#+M+l) comprises the new combined RAKE output introduced into the filter. 

As discussed above, the correlation matrix for the solution for the q th symbol may be 
represented by: 
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R z (q) = E{z q z»}, (58) 

where E{-} is the statistical expectation operation, which in this case is taken only over the 
random noise. Substituting Equation 56 into Equation 58 results in: 



R,(<7) = £ 

Repeating this with Equation 57 for q + 1 results in: 

R 2 (9 + i) = £ 



z(q-M)z*(q-M) z(q-M)a% 
a q z*(q-M) 



a a" 

11 



(59) 



a a" 



a z(q +M + 1) 



z*(tf + M + l)a^ z(tf + M + l)z*(? +M + 1) 

As shown in Equations 59 and 60, a submatrix is {a^a" } appears in the correlation 
matrix for symbol q and q + 1, indicating that it may be reused. Essentially, for each 
successive symbol sub-matrices ^{a^a" } may be shifted to partially form a new correlation 

1 0 matrix for the new symbol. Figures 9a and 9b illustrate the nature of this matrix shift. As 
shown in Figure 9a, 600 comprises row one 610, column one 620, and submatrix 
E^i q 2i" }630. Row one 610 and column one 620 border the leftmost and uppermost 

boundaries of submatrix ^ja^a" } 630. Further, Figure 9b illustrates that R^+i 602 comprises 

row M 612, column M 622, and submatrix £"{a 9 a" } 630. R^+i 602 may be generated from 

15 by shifting submatrix a q a q H 630 into the upper left-hand corner of R^ +i 602. As a result, 

only row M 612 and column M 622 is calculated to complete R^+i 602. Depending on the 
size of R matrix, the computational savings may be substantial. 

The MMSE solution allows for the special time varying structure of the data 
correlation matrix, R 2 . As noted above, the variation of the spreading waveform necessitates 

20 the calculation of a new weight vector for each symbol period. This variation of the 

spreading waveform manifests itself in the changing of h' 0 and R/ 1 . However, R z does not 
completely change from symbol to symbol, and therefore, a substantial portion of R z may be 



(60) 
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reused. As a result, only a portion of R z needs to be calculated from scratch for each new 
symbol. 

Further, in the multi-code scenario, the same MC-LEQ 500 with the same data 
contents may be used to solve for all K code channels at a given symbol instance. If the ML 
5 solution is used, the K u matrix is recalculated for each code channel because R w varies from 
code channel to code channel. However, with the MMSE solution, a single common data 
correlation matrix, R 2 may be calculated and used for all K solutions. 

The RAKE receiver embodiments discussed herein may be disposed in any wireless 
communication device, such as a base station 700 or mobile terminal 710, as shown in Figure 

10 1 0. As used herein, the term "mobile terminal" may include a cellular radiotelephone with or 
without a multi-line display; a Personal Communications System (PCS) terminal that may 
combine a cellular radiotelephone with data processing, facsimile and data communications 
capabilities; a personal data assistant (PDA) that can include a radiotelephone, pager, 
Internet/intranet access, Web browser, organizer, calendar, and/or a global positioning system 

1 5 (GPS) receiver; and a conventional laptop and/or palmtop receiver or other appliance that 

includes a radiotelephone transceiver. Mobile terminals may also be referred to as "pervasive 
computing" devices. 

Base station 700 includes antenna 702 for transmitting symbols to mobile terminal 
710. Objects, such as interfering object 720 may cause multiple echoes of the transmitted 

20 symbols to arrive at mobile terminal 710 at different times, as described above. The symbols 
received at antenna 712 of mobile terminal 710 are then processed in RAKE receiver 714. 
RAKE receiver 714 represents any of the exemplary embodiments described above. 
Similarly, mobile terminal 710 may transmit symbols along multiple paths to base station 
700. RAKE receiver 704 in base station 700 may then process the received symbols 

25 according to any of the exemplary embodiments described above. 
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In discussing exemplary embodiments of the present invention, it should be 
understood that one or more embodiments of the present invention comprise signal 
processing methods that may be implemented in hardware using discrete or integrated 
circuits, in software as stored program instructions, or in some combination thereof. More 
5 generally, one or more embodiments of the present invention may be embodied in hardware 
and/or software (including firmware, resident software, micro-code, etc.), including an 
application specific integrated circuit (ASIC). 

The present invention may, of course, be carried out in other ways than those 
specifically set forth herein without departing from essential characteristics of the invention. 
1 0 The present embodiments are to be considered in all respects as illustrative and not 
restrictive, and all changes coming within the meaning and equivalency range of the 
appended claims are intended to be embraced therein. 
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